Optical parametric oscillation with distributed feedback in cold atoms 
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There is currently a strong interest in mirrorless las- 
ing systems [1], in which the electromagnetic feed- 
back is provided either by disorder (multiple scat- 
tering in the gain medium) or by order (multiple 
Bragg reflection). These mechanisms correspond, re- 
spectively, to random lasers [2] and photonic crystal 
lasers [3]. The crossover regime between order and 
disorder, or correlated disorder, has also been inves- 
tigated with some success [4—6]. Here, we report one- 
dimensional photonic-crystal lasing (that is, distrib- 
uted feedback lasing [7, 8]) with a cold atom cloud 
that simultaneously provides both gain and feed- 
back. The atoms are trapped in a one-dimensional 
lattice, producing a density modulation that creates 
a strong Bragg reflection with a small angle of inci- 
dence. Pumping the atoms with auxiliary beams in- 
duces four-wave mixing, which provides parametric 
gain. The combination of both ingredients generates 
a mirrorless parametric oscillation with a conical out- 
put emission, the apex angle of which is tunable with 
the lattice periodicity. 

Among the possible systems that can be used to produce 
and study mirrorless lasers, cold atoms are interesting be- 
cause of their specific properties that differ from these of 
standard photonic materials. First, they are resonant point- 
like scatterers, producing extremely narrow spectral features 
(gain curves, scattering cross-section), which can provide flex- 
ibility [9] or give new effects. Second, their temperature is low 
enough to make Doppler broadening negligible in most situa- 
tions, but large enough to make them move substantially on 
a millisecond timescale, which makes disorder-configuration 
averaging or dynamic evolution from order to disorder very 
easy. Third, cold atoms are well isolated from the environ- 
ment, which makes them good candidates in the search for 
quantum effects. 

Conventional lasing has already been demonstrated when 
cold atoms are used as the gain medium [10, 11], as has radi- 
ation trapping due to multiple scattering [12], and efforts are 
under way to combine both factors to obtain random lasing 
[13,14]. In the opposite regime, a one-dimensional photonic 
bandgap (PBG), yielding efficient Bragg reflection of light, 
has recently been demonstrated in a cold, ordered atomic 
vapour [15]. In this Letter, we demonstrate optical paramet- 
ric oscillation (OPO) with distributed feedback (DFB) in cold 
atoms trapped in a one-dimensional optical lattice, by com- 
bining the PBG with four- wave mixing (FWM), which pro- 
vides the gain mechanism [10, 16-18]. 

We trapped cold 87 Rb atoms in a one-dimensional lattice 
of tunable wavelength Adi P . The trapping beam is retrore- 
flected to generate a potential of periodicity Adi P /2 (Fig. 1 
and Methods). Typically, N = 5 x 10 7 trapped atoms were 
distributed over a length L « 3 mm (~ 7, 700 atomic layers) 



at a temperature T « 100 /xK, leading to a root-mean-square 
(r.m.s.) transverse radius of the cloud a± « 60 /xm. Such an 
atomic pattern gives rise to a periodic modulation of the re- 
fractive index n and we have shown recently [15] that the 
very small modulation amplitude An ~ 10 -3 inherent to di- 
lute vapours (density p ~ 10 12 cm -3 ) could be balanced by 
the large number of layers, provided that the Bragg condition 
is fulfilled for a small angle of incidence and for a frequency 
slightly off the atomic resonance to avoid too much scattering 
losses. A Bragg reflection as efficient as 80% can be obtained. 

In the present experiment, we investigated the situation 
when gain was added to the system. Cold atoms can amplify 
light when pumped by auxiliary near-resonant beams, and 
several gain mechanisms have already been demonstrated (see 
ref. [10] and references therein). The combination of enough 
gain and multiple Bragg reflection should lead to DFB lasing. 

However, in this system, the stability of the feedback mech- 
anism is a critical issue. We trap the atoms using a lat- 
tice beam that is far red-detuned from the atomic transition 
(D2 line of 87 Rb, wavelength A = 780.24 nm m vacuum, 
linewidth r/27r = 6.07 MHz), so that the Bragg condition 
can only be fulfilled for a non-zero propagation angle 0, rel- 
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Figure 1. Schematics of the set-up. Cold atoms are trapped 
in the lattice formed by a retroreflected dipole trap. The pump 
beam is also retroreflected and has an incident angle of ~ 8°. 
Above threshold, the system emits light with an angle around 
the lattice: in any given plane including the lattice, four waves are 
coupled. Ei and E4 as well as E2 and E3 are coupled by the phase- 
conjugation process; E\ and E3 as well as E2 and E4 are coupled 
by the Bragg reflection. The emitted light is detected by avalanche 
photodiodes (APDs) and the beam cross section is observed by a 
charge-coupled device (CCD) camera. Additionally, a probe beam 
(incident power Po) can be used to measure the phase-conjugate 
reflectivity R pc = P pc /Po, where P pc is the reflected power. Inset: 
scheme of the four-photon transition corresponding to FWM. Up- 
ward (downward) arrows represent pump (probe and conjugate) 
photons. MOT, magneto-optical trap. 
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Figure 2. Transverse mode of the emitted beam, a, Cross-section of the beam for different Adi p . Each image is the result of 
an average of about 100 pictures of integration time 200 /is each, b, Measured emission angle 6 as a function of Adi p (red symbols). 
Errors bars indicate the uncertainty in fitting the radius on images and the angle calibration (see Methods). The dashed line is the Bragg 
condition (equation 1) with n = 1 and the solid line is a fit with n as a free parameter, c, Angular correlation function C(Acp) computed 
from the Adi p = 780.75 nm data. The solid line is the average over all images and the blue area repressents the statistical uncertainty. 



atively to the lattice axis, given by 

Adip cos = Xo/n (1) 

where n is the averaged refractive index of the atomic 
medium. Because of this angle, typically # ~ 2° for Adip = 
780.6 nm, the light beam is not perfectly reflected onto itself 
such that it transversely walks off and leaves the interaction 
volume after a certain number of Bragg reflections. 

Because this walk-off effect plays an important role, we 
must consider, for each plane containing the lattice, four 
different waves having a propagation angle relatively to 
the lattice axis (Fig. 1). In that case, a 'standard' gain am- 
plifies one wave (e.g. E4 in Fig. 1) while the Bragg reflec- 
tion couples it to another one (E2) with a walk-off. This 
problem can be overcome by producing gain with a phase- 
conjugation mechanism such as degenerate or nearly degener- 
ate FWM [10,16-18]. Then, each wave generates a backward 
phase- conjugated wave. Combined with the Bragg reflection, 
this leads to a global, walk-off free coupling between all four 
waves (Fig. 1), which favours an oscillatory behaviour. 

Experimentally, inducing FWM in a phase-conjugation 
configuration is done by simply retroreflecting a near- 
resonant, linearly-polarized pump beam (Fig. 1). Its typi- 
cal detuning from the F — 2 — >• F' — 3 closed transition is 
A = — 5r, it makes an incident angle with the lattice of ~ 8° 
and it is collimated with a waist w ~ 2.4 mm, thus ensuring a 
nearly- homogeneous pumping of the whole lattice. To avoid 
optical pumping into the dark hyperfine F=l ground state, 
a repumping laser is kept on all the time. In addition, a weak 
probe beam, phase-locked with the pump, can be used for 
pump-probe experiments or as a local oscillator. 

When the pump power P overcomes a certain threshold, 
namely Pth ~ 1 mW for A — oT (Fig. 3), we observe a 



strong, directional light emission that can be recorded either 
with avalanche photodiodes or with a charge-coupled device 
(CCD) camera (Fig. 1). This radiation is due to OPO with 
DFB in the cold atom sample, with FWM as the gain mecha- 
nism. This interpretation is supported by many observations, 
as discussed in the following. 

First, this radiation is obtained only with a retroreflected 
pump beam, whose alignment is critical, which is a strong in- 
dication that FWM is at work (with only one pump beam, we 
observe a strong Raman gain [10] in pump-probe experiments 
but no laser). Second, the polarization of the emitted radia- 
tion is linear and orthogonal to the one of the pump beam. 
This is also consistent with the properties of FWM, which is 
much more efficient for orthogonal pump and probe polariza- 
tions [19]. Third, we measured the frequency of the emitted 
light by a beat note experiment (see Supplementary Informa- 
tion). The emitted light is just a few kilohertz detuned from 
the pump frequency, which is consistent with nearly degener- 
ate FWM [20] and inconsistent with Raman gain, for which 
the amplification line was detuned by ~ 200 kHz with similar 
parameters. 

The role of Bragg reflection as the feedback mechanism 
is demonstrated by the beam shape. Emission occurs with 
angle from the lattice axis and, because of axial symme- 
try, forms two cone-shaped beams on each side of the lattice. 
From images of the beam cross-section (Fig. 2a), we extract 
the emission angle as a function of Adi P . A fit with ft as a free 
parameter is fully consistent with the Bragg condition (equa- 
tion 1) and gives n - 1 = (2.2 ± 0.5) x 10~ 4 (Fig. 2b), which 
is the expected order of magnitude given the atomic density 
(the refractive index depends also on the pump power). 

The intensity profile I ((f), where (p is the in-plane angle, 
exhibits strong shot-to-shot fluctuations. However, there is al- 
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Figure 3. Emitted laser power and maximum FWM re- 
flectivity as a function of the pump power. Red filled 
squares: laser intensity. Blue filled circles: FWM reflectivity below 
the threshold. The pump detuning is A = — 5I\ Blue open cir- 
cles: FWM reflectivity measured without lattice (to compensate 
the lattice-induced light-shift, we increased the pump detuning to 
A = — 7r). Inset: zoom-in for low pump power. The comparison 
between the reflectivity measured with and without the lattice 
shows that the lattice increases the reflectivity when approach- 
ing the threshold, which demonstrates a positive feedback effect. 
Above threshold, only the reflectivity without lattice can be meas- 
ured. 



ways a symmetry, that is I(ip) « I {if + 7r). This can be more 
precisely quantified by computing from the images the angu- 
lar correlation function C(Aip) = (I(ip)I(ip + Aip)) / (I ((f) 2 ) , 
which shows indeed a very strong correlation C(tt) ~ 0.96 
(Fig. 2c). This correlation comes from the couplings, namely 
the Bragg reflection and the phase conjugation, that exist 
in any given plane between the four directions of emission 
(Fig. 1). In contrast, the shot-to-shot fluctuations can be un- 
derstood by the random direction of the initial spontaneous 
emission event that triggers the laser oscillation and by the 
absence of coupling between waves in different planes. In ad- 
dition, we note that even after averaging, the intensity is not 
uniformly distributed along the ring (Fig. 2a), indicating that 
some effects break the axial symmetry (the incident angle of 
the pump beam and possibly some residual astigmatism in 
the lattice beam). 

Finally, we have measured the emitted intensity as a func- 
tion of the pump parameters. We observed that the OPO is 
not very sensitive to the precise value of the pump detuning 
A, provided it is negative (red-detuned) and not too close to 
resonance: we did not observe any significant change between 
A = — 4T and A = — 10F. Closer to resonance or with a posi- 
tive detuning, the pump-induced heating destroys the lattice, 
and the laser emission is inhibited or much weaker. The de- 
pendance on pump power is much more important. Above 
the threshold, the emitted power increases with pump power 
up to an optimum at P ~ 5 mW, above which a decrease is 
observed (Fig. 3). We estimate the maximum emitted power 
to be ~ 3 fjiW on each side. The decrease of the reflectiv- 
ity past an optimum pump intensity is a known behaviour 
of FWM [16]. In our case, we also suspect a detrimental me- 
chanical effect (heating or residual radiation pressure) which 
destroys the lattice. We indeed observed, in the temporal be- 



haviour of the emitted light, that radiation is sustained for a 
longer time for lower pump intensities, the longest duration 
being ~ 0.5 ms. 

In summary, it is interesting to note that the threshold 
condition requires ~ 50% phase- conjugate reflectivity meas- 
ured with a disordered sample (Fig. 3) and ~ 80% Bragg 
reflectivity measured with a passive sample [15], but when 
both ingredients are combined, these quantities are modified 
and we can no longer determine the amounts of feedback and 
gain, in contrast to standard lasers, in which the gain medium 
and cavity can be characterized independently of one another. 
This is also what makes our system original and interesting. 

To conclude, let us discuss the possible applications of this 
DFB OPO. First, we stress that the idea of combining FWM 
with a red-detuned DFB grating to generate conical beams 
(and possibly Bessel beams [23]) might also be applied in 
other systems, such as semiconductors [21,22]. In contrast, 
on-axis feedback could be obtained by using higher-order 
Bragg reflection, for exemple, with Adi P = 2Xo/n. Combining 
this idea with frequency up conversion schemes [24-27], one 
might be able to produce mirrorless oscillation at Ao with 
pumps at longer wavelengths. In principle, Ao can be any 
atomic transition, possibly in the ultraviolet range. Another 
possible application is the generation of 'twin beams', that 
is, pairs of beams with a relative intensity noise below 
the standard quantum limit. It is indeed well known that 
FWM induces quantum correlations and it has been shown 
that FWM oscillators above threshold could create twin 
beams [20]. We have already demonstrated strong angular 
classical correlations, and more work is needed to investigate, 
theoretically and experimentally, the question of possible 
quantum correlations in such systems. Here, the correlations 
should concern the four directions of emission in each plane. 
Finally, extension to a three-dimensional geometry, using the 
bandgap predicted in ref. [28], is an exciting open question. 

Methods 

Atomic sample preparation and experimental cycle. The 

atoms were trapped and cooled in a magneto-optical trap (MOT) 
loaded from a background vapour in 1 s. Stages of compression 
and molasses were used to increase the density and decrease the 
temperature, before loading the atoms in the dipole trap, which 
was generated by a home-made continuous-wave Tirsapphire laser, 
following the design of ref. [29]. The maximum available power is 
1.3 W and its wavelength Adi p was tunable in the range 770 820 nm. 
The beam was focussed with a waist (1/e 2 radius) u>di P = 220 /im 
at the MOT position (Rayleigh length ~ 20 cm). The lattice was 
then formed by retroreflecting the beam. After loading the optical 
lattice, the molasses beams were switched off and a waiting time of 
20 ms allowed the untrapped atoms to escape. Then, we could either 
characterize the trapped sample by absorption imaging, perform 
pump-probe spectroscopy to measure transmission, Bragg reflection 
and phase-conjugate reflection spectra, or shine only the pump beams 
to observe the emitted light. When the phase-conjugate reflectivity 
was measured without lattice (Fig. 3), the lattice was switched off 
with a mechanical shutter in 130 fis and the atoms freely expanded 
for 1 ms before measurements. This time of flight was large enough 
to completely smooth out the ordered structure and small enough 
to keep the optical thickness constant. All the stages following the 
initial MOT loading lasted only a few milliseconds, so that the total 
cycle duration was not much longer than 1 s. The repetition rate 
was thus ~ 1 Hz, which allowed quick averaging over many realizations. 

Detection tools. To measure the phase-conjugate reflectivity, we used 
APDs and obtained spectra by sweeping the probe frequency with 
an acousto-optic modulator in double-pass configuration. Because the 
probe power fluctuated during the sweep, the recorded reflected inten- 
sity was divided by a reference intensity recorded simultaneously with 
another APD illuminated by a part of the probe beam. The relative 
sensitivity of both APDs was calibrated with a 10% uncertainty due 
to thermal drifts. 

To observe the transverse mode of the OPO (cross-section images 
shown in Fig. 2), we used a CCD camera. The beam was first reflected 
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off a mirror, which was pierced to allow the lattice beam to go through 
(Fig. 1), and then collimated and refocussed 50 cm later on the camera. 
This allowed us to image the whole beam, and also to use several small 
black masks mounted on thin glass plates to get rid of stray reflections 
of the lattice beam on the vacuum chamber windows, which focused 
at intermediate distances. The angle calibration used the probe beam, 
which could also be imaged on the camera, as a reference. The incident 
angle of the probe beam was precisely determined using a series of 
Bragg reflection spectra off the passive lattice. When the spectra were 
symmetric, the probe angle fulfilled cos 9 = Ao/Adi p [15]. 

Received 7 July 2011; accepted 14 November 2011; published 
online 18 December 2011. 
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Supplementary information 

Beat note experiment 

To determine the frequency of the emitted radiation, respectively to the pump frequency, we use our probe beam 
as a local oscillator and perform a beat note experiment. 

Both beams are phase-locked together through an injection locking, which allows the resolution of very narrow 
spectral features. However, both beams pass through different double-pass acousto-optic modulators (AOMs) to 
allow us to sweep their frequencies. Standard AOM drivers suffer from thermal drifts and thus do not provide a 
kHz precision. As a consequence, for this beat-note experiment, we use two signal generators (Rohde & Schwarz 
SMR20) to drive the AOMs with a well-defined frequency difference of fo/2 — 150 kHz (because of the double-pass 
configuration, the frequency difference between both laser is fo = 300 kHz). 

We overlap the probe beam on one of our detection photodiodes and record the temporal trace of the emitted 
radiation beating with the probe. We generate the Fourier transform of this signal with a digital oscilloscope and 
we average this power spectrum S(f) over many realizations. We observe a peak around fo (inset of Fig. S4), whose 
precise position gives the relative frequency shift 5f between the pump beam and the emitted radiation. The precision 
of the measurement is limited by the duration of the OPO. 

Even though the most efficient configuration for four- wave mixing is when all waves are degenerate, we observe a 
non-zero frequency difference. This can be explained by the role of the refractive index, which provides the feedback, 
and whose value is not necessarily maximum at the pump frequency. The resulting working frequency of the OPO 
comes thus from a trade-off between the gain efficiency and the feedback efficiency. 

We measured the frequency difference Sf for different pump powers and observe a linear drift (Fig. S4), which is 
consistent with a broadening of pump-induced spectral structures. 
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Figure 4. Beat note experiment. Red dots: frequency difference Sf between the pump and the emitted radiation, as a 
function of the pump power. We observe a linear frequency drift (dashed line). Inset: example of measured spectrum. 
The error bars of the main panel correspond to the peak width (FWHM) of the beat-note spectrum. 
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